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Applying a variational Monte Carlo method to a two-dimensional t-J model, we study the 
nonmonotonic d x 2_ y 2-~w&ve superconductivity, observed by Raman scattering and ARPES ex- 
periments in the electron-doped cuprates. As a gap function in the trial state, we extend the 
d-wave form (ext. a!) so as to have its maxima located near the hot spots of the system. It 
is found that, in contrast to the hole-doped case, the ext. a! wave is always more stable than 
the simple d wave in the electron-doped case, and the magnetic correlation of the wave vector 
(n,Tv) as well as the pair correlation is enhanced. These results corroborate spin-correlation- 
mediated superconductivity in cuprates, recently argued from a FLEX calculation. In addition, 
we confirm that s- and p-wave symmetries are never stabilized even in the over-doped regime. 
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Introduction: Antiferromagnetic (AF) correlation is 
probably the primary origin to form Cooper pairs in the 
high-T c cuprates. 1,2 A recent neutron scattering experi- 
ment 3 in the electron-doped (n-type) cuprates discovered 
peaks at magnetic Bragg spots in both normal and su- 
perconducting (SC) phases, which fact indicates that the 
spin correlation of the AF wave vector Q = (tt, it) plays 
an important role in n-type cuprates. 

As in the hole-doped (p-type) cuprates, the pair- 
ing symmetry of the n-type ones was ascertained to 
be a d a ,2_j / 2-wave type, as far as the doping rate S is 
smaller than 0.15, by a scanning SQUID microscope 4 and 
angle-resolved photoelectron spectroscopy (ARPES). 5,6 
However, recent experiments by Raman scattering 7 and 
ARPES 8 have concluded that the pairing symmetry in 
the n-type cuprates is not the typical d x 2_ y 2 wave char- 
acterized by Ak oc cos k x — cosk y , but exhibits a non- 
monotonic behavior. Namely, the maximum of Ak is 
located midway between the Brillouin-zone boundary 
(tt, 0) and the zone diagonal (tt/2, it/2). Furthermore, for 



T > T c , pseudogap-like behavior 9 or an AF gap 10 arises -p Q _ TJ.(1 



results of Raman spectral functions and spin susceptibil- 
ity etc. with the experiments. Since the electron corre- 
lation is not weak even in the n-type cuprates, the re- 
sults of FLEX, which is basically a weak-coupling the- 
ory, should be checked by complementary studies from 
the strong-coupling and low-carrier-density sides. In this 
paper, we study this issue, applying a variational Monte 
Carlo (VMC) method 12 to a t- J-type model. So far, tak- 
ing this approach, various aspects of the cuprates have 
been elucidated on a strong-coupling footing. 13 

A secondary interest of this paper is possible variation 
of the pairing symmetry from d to s wave in the over- 
doped regime of n-type cuprates, observed by tunneling 
spectroscopy 14 and measurement of magnetic penetra- 
tion depth. 15 Although a BCS-level calculation 16 sup- 
ports these experiments, one has to confirm this issue 
using a less biased approach. 

Formulation: We consider a two-dimensional t-J 
model, H = U t +Uj with H t = - 2 (i j)a U 3 VG(c\ a c ja + 
R.c.)V G and Hj = JY,(a){Si • Si' 



at this midway locus, in contrast to the p-type cuprates, 
in which the gap maximum and pseudogap behavior take 
place around (ir,0). 

These results can be explained by the fact that the 
loci of the gap maximum in the p and n types roughly 
coincide with their respective hot spots — the intersec- 
tion of the Fermi surface and the magnetic Brillouin zone 
boundary [See Fig. 1(a)]. This observation affords confir- 
matory evidence that the superconductivity in cuprates 
is induced by the spin correlation (or fluctuation) of the 
wave vector Q, which connects the hot spots. 

Theoretically, Yoshimura and Hirashima 11 recently 
treated this issue, applying a fluctuation exchange ap- 
proximation (FLEX) to the Hubbard model. They ob- 
tained a non-monotonic behavior of Ak and consistent 



n n n 3l 



iHj/4:), where 
), the value of J/t is fixed at 0.3, 



and Uj = t,t',t" or 0, according as the site i is a first-, 
second-, third-nearest neighbor or farther site of the site 
j, respectively. An electron-doped (more-than-half-filled) 
system can be treated with a t-J model as a less-than- 
half-filled case, by applying a particle-hole transforma- 
tion ct CT -> exp^Q-r,)^ with t'(t") -> -t'{-t"). Thus, 
we put t'/t > (< 0) and t"/t < (> 0) for n- (p-)type 
cuprates. Actually, we adopt typical values 17 of t'/t and 
t" /t given in the caption of Fig. 1. 

To this model, we apply a VMC method, which accu- 
rately treats the local correlation Vq. As a variational 
function for a SC state, a simple Gutzwiller-type wave 
function is used, 



|*sc) = 7>g|$bcs) = V G 
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Fig. 1. (a) Hot spots are compared between a hole-doped (t' /t = 
—0.1, t"/t = 0.1, open circle) and an electron-doped (t'/t = 
—0.16, t" ' /t = 0.2, solid circle) cases for the optimal doping (8 = 
0.15). The bare Fermi surface for the former (latter) parameter 
set is denoted by a thin-dashed (solid) line, and the magnetic 
Brillouin zone boundary by a bold dotted line. The AF nesting 
vector Q links two hot spots, (b) Pairing potential A k /t of the 
ext.d wave for the optimized parameter values (see text) of the 
system shown in (c). (c) The fe-points on the Fermi surface for 
an underdoped density on a 14 X 14 lattice are shown with solid 
circles. Regions near the hot spots are indicated by shadows, and 
the magnetic Brillouin zone boundary by a dashed line. 



with 



fk = 



Ak 



e k - n + V( £ k"M) 2 + A k ' 



(2) 



Here, $bcs is the BCS wave function of a fixed particle 
number N c , £k = — 2i(cos k x + cos k y ) — At 1 cos k x cos k y — 
2t"(cos2fc a; +cos2fcj / ), and the parameter fi is substituted 
by the chemical potential of the non-interacting case. It 
is already known that this type of wave function works 
well for i-J-type models. 18 ' 19 Anisotropy of the pairing 
potential is introduced into Ak, as A k = Ar^cosfcc — 
cos k y ) for the simple d wave. We extend it so that it 
may have large amplitude near the hot spots, 
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and call it the extended-d (ext.d) wave. In A k xt d , we 
add to A k eight Gaussian peaks, whose height, position 
and width in two directions are adjusted by parameter 
a,b,c and d, respectively. Although Ak should be deter- 
mined variationally, for simplicity we fix the parameters 
for the Gaussian, except for trial calculations, 20 at rough 
optimal values of the case 5 — 0.082 (L — 14), namely, 
a = 3, 6 = 2, c = 1 and d = 1; Ak thereof is shown 
in Fig. 1(b). Note that the loci of the maximal Ak are 
situated very closely to the hot spots, shown in Fig. 1(c). 
Thus, Ao (amplitude of Ak) becomes a sole parameter 
to be optimized. Unlike the BCS theory, Ak here is not 
entirely equivalent to the SC gap, especially in A , but 
the symmetry of Ak faithfully reflects the SC gap. 

We have followed a conventional VMC scheme, 12 and 
collected I0 5 -I0 6 samples, which suppress the statistical 
errors in energy at ~ 10~ 4 t. The system used has L x L 
sites (L = 12, 14 and 16) with the periodic-antiperiodic 
boundary conditions. The particle density 5 is chosen so 
as to satisfy the closed-shell condition. 
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Fig. 2. Total energies versus a variational parameter Ao for var- 
ious pairing symmetries for an electron-doped system with S = 
0.163. The arrows on the d and ext.d waves denote the optimal 
values. The value of the normal state (Ao = 0) is shown with a 
dashed line. The system is of 14 X 14, and J/\t\ = 0.3. 



Results: Before going to the d-type pairings, we check 
the stability of the singlet s- (Ak = Ao), ext. s-[= 
A (cosfca; + cosky)] and triplet p-wave (= A n smk x ) 
symmetries. In Fig. 2, the variational energies E tot of 
the symmetries we treat are compared in the overdoped 
regime. Here, E tot 's of the s, ext.s and p waves mono- 
tonically increase, as A increases, and are destabilized 
with respect to the normal state. We have also confirmed 
that the same behavior persists to a large doping rate, 
S = 0.245, where SC states are no longer stabilized. Since 
our treatment is little biased in comparing pairing sym- 
metries, we are confident that the p and s-type waves 
are not realized even in the overdoped regime; thus a 
pairing-symmetry transition is unlikely to arise. 21 

Now, we turn to the d-type symmetries. By contrast, 
the variational energies of the d- and ext.ci-wave SC 
states plotted in Fig. 2 have minima at finite values of 
Ao, as indicated by arrows. The decrease in E to t of the 
d x 2_ y 2 wave has been well-known for the plain t-J model 
(t 1 = t" = 0) since the early stage. 25 It should be noted 
here that the ext.d wave has an appreciably lower energy 
than the simple d wave even for such a large value of 5. 
We have optimized E to t similarly for various values of 5, 
and depict in Fig. 3(a) the difference of E tot between the 
d and ext.d waves, AE = E? ot - £ t °* t t d . 26 In electron- 
doped cases, the ext.d wave is always more stable than 
the d wave. The large value of AE near half filling proba- 
bly stems from the fact that the Gaussian peaks in A k xt d 
are close to the hot spots, as well as that the energy scale 
in the condensation energy E cond (= £ t n ° t rmal - Ef£) be- 
comes large for S — > 0. AE vanishes at 6 = 0.222, where 
the hot spots still survive, but E cond vanishes for both 
waves. Conversely, in the hole-doped cases [open symbols 
in Fig. 3(a)], the simple d wave is more stable than the 
ext.d wave except in the vicinity of half filling. In the 
optimal- and overdoped regime for S < 0, the hot spots 
sit near (it, 0) and equivalent points, whereas near half 
filling the hot spots are still away from (tt, 0) point and 
rather closer to the Gaussian positions we set [Fig. 1(b)]. 
Thus, the results in E to t definitely indicate that the SC 
state becomes stable when the maximum of the gap Ak 
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Fig. 3. (a) Difference in the optimized energy between the ext.d- 
and simple d-wave states, AE = Ef ot — Ef^- d , as a function 
of doping rate 5. Since the parameters (t' and t") are different 
between p and n types, the values at 8 = do not coincide, (b) 
Comparison of the pair correlation function between the ext.d 
and simple d waves as a function of <5. In both panels, hole-doped 
systems are indicated by negative values of S. 



is located near the hot spots of the system. 

To reinforce the above argument, we consider the d- 
wave SC correlation function of nearest-neighbor pairs, 

1 



j T,T 



£ (-l) 1 -^')(A T (j)A T ,(j + r)), 



with At(j) = (cj T c] +ri +c]+rT c jl)- Sincc Pd ( r ) ra P idl Y 
decays with |r| and is almost constant for |r| > 3, the 
average for |r| > 3, P^ ve , gives an adequate estimate of 
the long-distance value. In Fig. 3(b), P| vc is plotted ver- 
sus carrier density. The remarkable asymmetry between 
the p and n types can be attributed to the difference in 
DOS at the Fermi surface, 27 particularly at the hot spots. 
Note that the ext.d wave always exceeds the d wave for 
the n type, whereas the relation is inverse in the over- 
doped regime for the p type. This tendency of Pd( r ) cor- 
responds well with that of AE [Fig. 3(a)]. Incidentally, in 
n-type systems, the SC correlation of long-distance pairs 
is probably enhanced, as will be discussed later. 

Next, to identify the origin of the energy gain in E tot , 
we compare the energy components, namely hopping en- 
ergy E t — (Tit) and exchange energy Ej = (Hj), among 
the d-wave, ext. d-wave and normal states. In Table I, we 
list the raw values of the d-wave, Ef and Ej, and the dif- 
ferences between them and those of the other two states, 
namely, AE?(= Ef - fi"™ 31 ), AE t (= Ef - E° xt - d ), 
etc. for four kinds of doping. As compared with the nor- 
mal state, the d- (also ext.d-) wave SC state is stabi- 
lized by the noticeable decrease in Ej; conversely, E t is 
more or less increases. Such magnetic origin of supercon- 
ductivity is characteristic of t-J-type models, 18 ' 28 and 



Table I. Comparison of the energy components, Et and Ej, for 
four carrier densities with J/t = 0.3 and t being the unit. The 
energy components of the simple d wave, Ef and Ej, are entered 
in the first two lines. Entered in the middle and lower two lines 
are the differences of Et and Ej between the d-wave and the 
normal state and between the d and ext.d waves, respectively. 
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Fig. 4. Comparison of the momentum distribution function be- 
tween the ext.d and pure d waves along the path indicated by an 
arrow in each inset, for the electron-doped system with <5 = 0.122. 
In each inset, solid dots and a broken line denote the outmost 
occupied fe-points (the Fermi surface) due to and the AF Bril- 
louin zone boundary, respectively. The system is of 14 X 14. 



closely related to the kinetic-energy-driven SC mecha- 
nism in the strong-correlation regime of the Hubbard 
model. 29 In comparison of the d and ext.d waves, Ej xt d 
is lower whenever the ext.d wave has a lower total energy 
(5 <; — 0.06), whereas the ext.d wave always possesses 
somewhat higher E t . Thus, the gap structure of the ext.d 
wave has further advantages to gain magnetic energy for 
n-type cases over the simple d wave. 

Now, we consider the momentum distribution func- 
tion, n(k) = 1/2 X^cr( c kcr c k(T), to actually observe the 
gap behavior in the momentum space — a milder slope at 
quasi-fcp. In Figs. 4(a) and (b), we depict n(k) for an 
electron-doped case 30 along two paths, namely OP [in 

(a) ], which goes away from the hot spots, and QR [in 

(b) ], which penetrates the hot-spot area. In (a), n(k)'s 
for the d and ext.d waves exhibit similar behavior, and 
the d wave seems slightly mild. On the other hand, in (b) 
the ext.d wave makes an obviously milder curve around 
the hot spot. Thus, the gap behavior around the hot spot 
is enhanced in the ext.d-wave state. Incidentally, in the 
node-of-gap direction, (0,0)-(7r, 7r), the difference between 
the d and ext.d waves is very small, and a clear Fermi 
surface [discontinuity in n(k)] can be seen (not shown). 

Finally, we consider the spin correlation func- 
tion. In Fig. 5(a), the spin structure factor, 5(q) = 
WsEije 4q - (i - j) S(iJ) with S(i,j) = (S?Sp, of the d 
and ext.d waves is plotted for several densities of elec- 
tron doping. Both waves have the maximal amplitude at 
(it, 7r) for all the electron densities, which is consistent 
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Fig. 5. (a) Comparison of the spin structure factor between the 
ext.d and pure d waves along the path (0, 0) — > (tt, 0) — * (n, n) —* 
(0, 0) for four densities of electron doping. Plotted in the inset is 
the difference of <S(q) between the ext.d and pure d x 2_ y 2 waves, 
A5(q) = 5 ,cxt d (q) — S d (q). (b) Ratio of real-space spin correla- 
tion function between the d and ext.d waves, S ext - d (i, j)/S d (i, j), 
is plotted for several values of i — j| as a function of 8. The data 
shown are the average among the four (or eight) equidistant sites. 
The system is of 14 X 14. 

with the neutron experiment. 3 As shown in the inset of 
Fig. 5(a), the difference between the d and ext.ii waves 
is almost restricted to the vicinity of (it, n), where S'(q) 
of the ext.d wave is sizably enhanced. This enhanced AF 
correlation naturally leads to the energy gain in Ej, men- 
tioned above. Shown in Fig. 5(b) is the ratio of the real- 
space spin correlation function for the ext.d wave to that 
for the simple d wave. Although the d and ext.d waves 
exhibit almost the same values for the nearest-neighbor 
sites (|i— j| = 1), 5(i,j) of the ext.d wave for farther dis- 
tances considerably increases, especially in the x (or y) 
direction. This is because the harmonics of the d wave, 
cosnk x — cosnk y (n > 2), give substantial contribution to 
^ext.d li Unlike the pure d wave (n = 1), they elongate 
the coherence length (and magnetic correlation length); 
thereby the correlation strength in the electron-doped 
systems becomes effectively weaker. 31 Such tendency has 
been actually observed by various experiments. 32 

Summary: Using a variational Monte Carlo method for 
a t-J model, we have studied a nonmonotonic (ext.) d- 
wave superconducting state, in which the amplitude of 
the gap parameter At is intentionally enhanced around 
the hot spots, so as to agree with recent experiments of 
Raman scattering and ARPES. This ext.ci-wave state has 
an appreciably lower energy than the simple d-wave state 
for all the densities of electron doping (also very low hole 
doping). This stabilization of the ext.d wave is caused 
by the gain in magnetic exchange energy, accompanied 
by a marked increase in the spin correlation of the wave 
vector Q = (tt, tt). In addition, we have shown that the s- 
type and p waves are unlikely to take place in the high-T c 
regime of 5. Our results using a strong-coupling approach 
basically agree with the recent FLEX study; 11 thereby, 
it is ensured that the AF spin correlation plays a crucial 
role for the high-T c superconductivity. 

We believe that the essence of the nonmonotonic d- 
wave gap is grasped in this work, although we have sim- 
plified A£ xt - d by both a naive assumption of a Gaussian 
form and fixing the parameters controlling the Gaussian. 
We should address quantitative refinement as well as re- 



lated issues. For instance, (1) simultaneous optimization 
of all the variational parameters, including the renormal- 
ization of the quasi-Fermi surface, enables us to follow 
the continuous evolution of Ak versus 5. (2) The rela- 
tion to the AF ordered state is very important. 
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